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POWER SERIES SOLUTIONS OF 


THE THIELE -BUREAU REGULARIZED PLANAR 
RESTRICTED THREE-BODY PROBLEM 


Regularizing transformations which remove the singularities at the attracting 
masses from the equations of motion of the planar restricted problem of three - 
bodies are well known. They have been developed by Thiele, Burrau, Levi-Civita, 
Birkhoff, Lemaitre, etc./ 3 ) to study the solution of the equations of motion at a 
collision and to facilitate numerical integration of collision orbits. 

The planar restricted problem of three bodies constrains two primary 
masses to move in circular orbits about one another while a third body of infini- 
tesimal mass moves in their plane of motion and exerts no force on the primaries . 
Regularization calls for the elimination of one or both of the two singularities r j -1 
and r J 1 , where r l and r 2 are the distances between the two primaries and the 
third body. By simultaneous, or global regularization, we mean the process 
which removes both singularities with one transformation, whereas by local 
regularization we mean the removal of only one singularity by a single transfor- 
mation. 

Global regularizations generally result in equations of motion which are 
considerably more complicated than either the original equations or the equations 
obtained when a local transformation such as that of Levi-Civita is used. To 
effect a solution by power series, the need for higher-order derivatives exists. 
Szebehely ^has proposed a general method based on differentiating the equation 
of motion which is straightforward but somewhat cumbersome. It is particularly 
inadequate when an optimal power series integration technique, such as that of 
Van Flandern^ 5 ) is desired. With this method, the integration step size at each 
step is optimized with respect to computing time so that a large number of high- 
order derivatives may be needed. 

The global transformation proposed by Thiele (1895) and Burrau (1906) has 
advantages over the other global regularizations since the equations of motion in 
this system contain only derivatives and transcendental functions of the dependent 
variables. Since the trigonometric and hyperbolic functions in the equations 
themselves satisfy simple differential equations, the Thiele-Burrau regulariza- 
tion is particularly amenable to the calculation of the series coefficients by re- 
current power series techniques. 


THIELE -BURRAU REGULARIZATION 

The equations of motion in a coordinate system with the origin at the mid 
point between the primaries and rotating with the mean motion are (4) 
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and where /a is the ratio of the mass of the smaller primary to the sum of the 
two primary masses. Here and r 2 represent the distances between the pri- 
maries and the infinitesimal third body, and assuming that the more massive 
primary is to the right of the origin we obtain 


Tj - \z~%\ (?) 

r 2 " I z + %\ . 

The first integral of the equations of motion is 

J = 20> - \z\ 2 (3) 

where J is the Jacobi constant of motion. 


To regularize the equations of motion, transformations 
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are introduced into Equation (1), where w = fi + i v and where the function f(w) 
is required to be such that the singularities of the equations of motion are re- 
moved. The result of this procedure is< 4 > 
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where 


w‘ = dw/dr , 
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The global regularizing transformation of Thiele-Burrau is defined by 
Equation (4) with 

z = f(w) = ^ cos w (8) 

We then find immediately from Equations (2), (6), and (7) 

■j(cosh v - cos u) (9) 

(cosh v + cos u) (10) 

= ^(cosh 2 v - cos 2 u) (11) 


df 
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and it follows that 

u"~2r 1 r 2 v' = ^ sin u [8(2 fi - 1)- 4 J cos u 

+ (2 fi - 1) coshv(3cos 2 u - cosh 2 v) + 2 cos 3 u] (16) 

v M +2r 1 r 2 u' = sirih v [ 8 - 4 J cosh v 

+ (2 fi - 1) cos u (3 cosh 2 v - cos 2 ) + 2 cosh 3 v] (17) 
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which are the globally regularized equations of motion for the infinitesimal third 
body under the transformation of Thiele -Burrau. The dimensionless coordinates 
x, y and time t may be recovered by the formulas 

x “ \ cos u cosh v (18) 

y s ~ ‘I' sin u sinh v (19) 


t 



(cosh 2 v - cos 2 u) dry. 


(20) 


while the dimensionless velocities may be obtained by applying Equation (5) to 
Equations (18) and (19). The result is 


(v 1 cos u sinh v - u 1 sin u cosh v) 
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sin u cosh v + u 1 cos u sinh v) 
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RECURRENT POWER SERIES SOLUTION 


The solution of Equations (16) and (17) is assumed to be represented by the 
power series 


u ~ 



(23) 



( 2 *) 


where the coefficients are to be determined recursively. 
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Then Equations (16) and (17) may be written as 

2u" + (p - q) v' - (2fi - 1) a + a i 
2v" - (p - q) u* “ c + eg . 


Further, from Equations (2 5) -(2 8) we have the relations 


bu' 
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Substituting the series expansions into the auxiliary relations (35)~(38), (31), (32) 
and into the equations of motion, (33) and (34) and equating coefficients of equal 
powers of r we arrive at the recurrence relations 

i 

1 a i = J2 J u j bi-j (39) 
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i t> A = - T j ^ a^j (40) 
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for i > 1, and 
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for i > 0. It is to be noted that for optimal computing advantage may be taken 
of the symmetry in Equations (43) and (44) to nearly halve the arithmetic opera- 
tions involved in these recursions. 

Initial conditions of the problem yield u 0 , v 0 , ti* and v r Then 

a Q - sin u 0 

b 0 = cos u 0 
c„ = sinh v 0 
d 0 = cosh v Q . 

With these values at hand, proceeding through the above algorithm will determine 
all higher-order power series coefficients in terms of the preceding coefficients, 
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so that the solution may be extended optimally by analytic continuation using a 
variable step size and a variable number of terms in the power series expansion 
for each integration step. 

This power series method for the solution of the Thiele -Burrau Regularized 
Planar Restricted Three-Body Problem was programmed in double precision on 
an IBM 360-91 computer using the optimal integration procedures outlined in 
Reference (1) and compared with the non-regularized recursive power series 
method of the same paper. As an indication of applicability, the periodic orbit 
of Figure 1 (2) in the Earth -Moon system (n> ~ 1/82.45) with period 

P = 6. 1921693313196 


and initial conditions 


x 


0 


- . 7121285627653111 



y 0 - 


0, 


x 0 = 0. 

y 0 = 1.04935750983033 

reveals a 33% decrease in computing time as well as two more digits of accuracy- 
in the Jacobi constant after one complete orbital integration. 

Thus, although the full power of this integration technique will be realized 
with collisions or near collisions involving both primaries, the above example 
demonstrates remarkable results in less extreme problems. 


p 
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